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Abstract 

The generalized open XXZ model at q root of unity is considered. We review how associated 
models, such as the q harmonic oscillator, and the lattice sine-Gordon and Liouville models are 
obtained. Explicit expressions of the local Hamiltonian of the spin | XXZ spin chain coupled to 
dynamical degrees of freedom at the one end of the chain are provided. Furthermore, the boundary 
non-local charges are given for the lattice sine Gordon model and the q harmonic oscillator with open 
boundaries. We then identify the spectrum and the corresponding Bethe states, of the XXZ and 
the q harmonic oscillator in the cyclic representation with special non diagonal boundary conditions. 
Moreover, the spectrum and Bethe states of the lattice versions of the sine-Gordon and Liouville 
models with open diagonal boundaries is examined. The role of the conserved quantities (boundary 
non-local charges) in the derivation of the spectrum is also discussed. 
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1 Introduction 



The ultimate goal when studying a physical system is the derivation of the corresponding observables. 
In the framework of integrable spin chains the first step towards such an aim is the diagonalization 
of the corresponding transfer matrix by means of the quantum inverse scattering method, introduced 
by the St. Petersburg school [UI21E]- The diagonalization process relies primarily on the existence of 
certain exchange relations emerging from the specific algebra that rules the model, defined by either 
the Yang-Baxter [3j HI [5] or the reflection equation [61 [7J. Prom the physical viewpoint Yang-Baxter 
and reflection equations describe the factorization of multiparticle scattering in the whole and half line 
respectively, a unique feature displayed by 2D integrable systems. Due to analyticity requirements 
imposed upon the spectrum certain constraints, known as Bethe ansatz equations [21 [8] arise, whose 
exact form depends explicitly on the choice of representation. The significance of Bethe ansatz equations 
rests on the fact that their solutions yield all the physically relevant quantities such as exact S matrices 
[21 El [10], thermodynamic properties [111 [T2l [T3] and correlation functions [31 [T3]. In addition it was 
recently realized [15] that in a particular limit string theories may be mapped to TV = 4 super Yang- Mills 
gauge theories, which in turn can be associated to special examples of integrable spin chains. Thus one 
may implement the powerful Bethe ansatz techniques [21 [8] , in order to derive non-perturbatively the 
physically relevant quantities. This remarkable realization put integrable spin chains and Bethe ansatz 
into the fore of recent advances in string theory. 

The main aim of the present article is the investigation of the spectrum of a generalized XXZ model, 
and associated models, with special open boundary conditions that preserve integrability. Historically 
the spin ^ XXZ chain with diagonal boundaries was investigated in [7J[T6l[T7], whereas the corresponding 
spin chain with non diagonal boundaries was just recently solved for q root of unity [18], and for generic 
values of q in [19]. The spectrum and Bethe ansatz for the spin s XXZ model with non-diagonal 
boundaries were also derived in |20[ 121] . In [22] the spectral equivalence between the Hamiltonians of 
the XXZ model with a non diagonal boundary and a novel open model (the asymmetric twin chain) 
with an obvious reference state is shown via the representation theory of boundary Temperley-Lieb 
algebras. Furthermore, the spectral equivalence of XXZ type Hamiltonians with diagonal and non- 
diagonal boundaries was investigated in [23], again in the context of boundary Temperley-Lieb algebras. 
Equivalence of the spectra of spin chains with different boundaries was also established in [21] by 
introducing appropriate defects. In the present study we shall be mainly interested in the case where 
the parameter q is a root of unity, and we shall restrict our attention in finding the spectrum of the 
general XXZ model and the q harmonic oscillator with non diagonal boundaries for the so called cyclic 
representation [25], which has not been treated so far. It will be instructive to provide a special example 
of the Hamiltonian of the spin ^ XXZ model with special integrable conditions, i.e. 

1 N ~ 1 N 
H = -- ^2 (°f°f+i + °f°f+i + cosh(^) ofof+i) - — cosh(z/i) + M N + ci erf + c 2 af + c 3 erf (1.1) 

i=l 

a x,y ' z , (<r ± ) are the 2x2 Pauli matrices. The constants Cj are dictated by integrability, and A4 is a 2 x 2 
matrix with entries being either scalars (c-number K matrix applied at the right boundary) or elements 
of a particular algebra, e.g. Ug{sh), the q harmonic oscillator algebra etc., that is a dynamical boundary 
coupled to the right end of the spin chain, that is a typical quantum impurity problem. In general such 
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type of Hamiltonians turn out to be of particular significance in condensed matter physics, describing 
for instance quantum impurities (see [131 EH 121! , [26]-[29] and references therein) or being related to 
the Azbel-Hofstadter Hamiltonian \26 \ \30 \ \31 \ I32j. Note that we shall also examine the spectrum and 
Bethe ansatz equations of the open lattice sine Gordon and Liouville models, but only with diagonal 
boundaries. 

Let us outline the content of each section of the present article. Sections 2 and 3 serve basically as a 
general review. More precisely, in section 2 we introduce the generalized XXZ model, and we also review 
how associated models such as the lattice sine-Gordon, Liouville models and the q harmonic oscillator 
arise via certain limiting processes. The underlying algebraic description is also provided and useful 
realizations are presented. In fact, one of the the main objectives of this article is to offer a unifying 
framework for the study of a whole class of open integrable lattice models emerging from the generalized 
XXZ model. However to identify the spectrum of the relevant models one has to fix a representation 
and this is done in sections 6 and 7. In section 3 section well known c-number K matrices, solutions 
of the reflection equation are reviewed, and the algebraic transfer matrix of an open spin chain is also 
recalled. 

In section 4 a special example of a XXZ Hamiltonian (jl.ip . with spin ^ in the bulk coupled to 
some dynamical system at the one end is considered. More precisely, explicit novel expressions of the 
dynamical boundary term M of the Hamiltonian (jl.ip are provided. Note that these expressions are 
generic (algebraic) at this stage, i.e. they are independent of the choice of representation for both the 
Mq(sh) an d the q oscillator algebras, which are our primary interest in the present study. In fact such 
a Hamiltonian may be regarded as a special example of the more general case described in section 
3. To deal with the corresponding spectrum naturally one has to choose a particular representation 
associated to the boundary. This is done however in section 6, where we restrict our attention to the 
cyclic representation for both XXZ and the q harmonic oscillator. In section 5 the boundary non local 
charges for the generic XXZ chain are reviewed, and novel expressions of the non-local charges for the 
lattice sine-Gordon model and q harmonic oscillator are derived. Their significance and their relation 
to the spectrum becomes clear in section 6. 

In section 6 the spectrum of the XXZ model and the q harmonic oscillator in the cyclic representation 
|25j . for special non-diagonal boundaries is investigated. Generalizing the formulation of [19\ 133] we are 
able to identify a reference state, by means of suitable local gauge transformations known also as 
Darboux matrices [M], upon which eigenstates are built. The derivation of the reference states, which 
is a new result, is achieved by solving sets of recursion relations. Within the algebraic Bethe ansatz 
framework novel expressions of the spectrum and Bethe ansatz equations of the relevant models are 
found. In this context the Hamiltonian presented in section 4, with the boundary terms associated to 
the cyclic representation, (see also (|l.ip ) is also treated as a special case. Another intriguing new result 
is the connection between the spectrum of the open transfer matrix and the spectrum of the conserved 
boundary non local charges, which is presented in section 6.3. Although this is a generic result, which 
hold for any representation we provide particular examples for the XXZ model and the q harmonic 
oscillator in the cyclic representation, given that these models are our primary concern. In Appendix 
C we present such a conserved quantity in the spin s representation (locally) as a tridiagonal (Jacobi) 
matrix |35j . and in the cyclic representation ofU q (sl2), and we attempt its diagonalization. In section 
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7 we give a flavor on the lattice sine-Gordon and Liouville theories with diagonal boundary conditions. 
This is the first time to our knowledge that these models with open boundary conditions are considered. 
Finally, in the conclusion section we summarize the results of the present article. 



2 The XXZ and associated models 



In this section the XXZ model is introduced, and various closely related models, which arise from it 
naturally are reviewed. The XXZ spin ^ R matrix, associated to the fundamental representation of 
Mq(sh)i acting on (C 2 )® 2 . In general, the R matrix is a quantity proportional to the physical S matrix, 
acts on V® 2 , and is a solution of the Yang-Baxter equation [3l [4] 

-Ri2(Ai — A2) i?i3(Ai) ^23(^2) = -^23(^2) #13(^1) Ru{^i — A2), (2.1) 

R<g>I, #23 = I <g) R. The spin \ XXZ R matrix in 



which acts on V , and as customary R\2 
particular is given by 



R(X) 



sinh/i(A + § + 
a + sinh(i/i) 



a smh(i/i) 
sinh / u(A + I - 



(2.2) 



The R matrix is written in the so called principle gradation, it can be however expressed in the homo- 
geneous gradation by means of a simple gauge transformation 



R { $(\) = Vi(-A) R ( f 2 \\) Vi(A), V(A) = diag(l, e» x ). 



? (p) 



(2.3) 



One may now derive a more general object L(X) E End(C 2 ) (g) U q {sl2), where the quantum algebra 
Mq(sh) is defined by the fundamental algebraic relation (see also Appendix A) 



Rabi^i — A2) L an {\i) Lb n {\2) — Lb n {\2) L an (\i) R a b(M — ^2J- 
A simple solution of the latter equation, which we shall use hereafter, is (index free notation) 



L(A) 



e^A-e-^D (q - g~ x )B 
(g-g -1 )C e^ A D-e^ A A 



(2.4) 



(2.5) 



A, B, C, D satisfy the defining relations of U q (sl2) (q = e l/i ), namely 
A D = D A = I, A C = gC A, A B = g~ x B A, 



C, B 



A 2 -D 2 



(2.6) 



A well known representation of U q {sl2) is the spin s representation, which is n = 2s + 1 dimensional, 
and may be expressed in terms of n x n matrices as 



n-1 



C — Ck efc k+i, 



k=l k=l 

where we define the matrix elements (eij)ki = $ik $ji and 

n + 1 



n-1 

B = Ck efc+i k 
k=i 



(2.7) 



fc, 



,-fc 



[k] q [n- k] q , [k]^ 



(2. 
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The generators A, B, C, D may be also expressed in terms of the Heisenberg-Weyl group generators 

X, Y 



XY = qYX (2.9) 



as 



A = D~ 1 =X, B= riq^'X' 1 -<fX)Y~\ C = -(q~ s X - q'X^Y. (2.10) 

q-q- 1 q-q- 1 

Such a realization may be also thought of as a g-deformation of the sfa algebra, when the corresponding 
generators are expressed in terms of differential operators, i.e. 

S z = y^~, S-=y- 1 (y-^ + s), S+ = y(y-f-a). (2.11) 
dy dy dy 

It is worth pointing out that the latter expressions for s = 0, correspond to a sl(2, R) realization, which 
has been employed in the context of high energy QCD. More precisely, it was shown in [36J that the 
high energy asymptotic behavior of the hadron-hadron scattering amplitude in QCD is described by the 
non-compact XXX Heisenberg chain with s = |37} [38] . 

In the present article the XXZ model will be examined for the special case where q is root of unity, 
i.e. q p = 1, q = e ZfJj , fx = where k, p integers. In this case the algebra admits a p dimensional 
representation, known as the cyclic representation [25]. More specifically, in addition to (|2.9j) one more 
restriction is applied so one may obtain a representation with no highest (lowest) weight 

X p = Y p = 1 

then the generators X, Y may be expressed as p dimensional matrices 

v p-i 
X = ^q- k e kk , Y = J2 e k k+i + epi. (2.12) 

k=l k=l 

This realization is of particular interest and it has been extensively used for instance in the problem 
of Bloch electrons in a magnetic field, described by the Azbel-Hofstadter Hamiltonian |30[ [31] 132] . In 
addition as argued in [39] the quantum Hall effect states may be also seen as cyclic representation of 
U q (sl 2 ). 

The XXZ model can be regarded as a universal model, since it gives rise to a whole class of associated 
integrable models. In what follows we shall briefly review how the lattice sine-Gordon and Liouville 
models are obtained in a natural way from the XXZ L matrix. Also the q harmonic realization will 
be obtained again from the generalized XXZ form. The generators X and Y may be associated with 
an infinite dimensional representation in terms of some lattice 'fields'. Then one may identify the L 
operator of the lattice sine-Gordon model [30]. More specifically set Z n = Y n X n , it is then clear that 
X n Z n = q7L n X n , where notice that for convenience we reintroduce the index n. Also set 



X n ► e ^ n , Wj n > e^ n 



i^5 mn (2.13) 
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q s ~2 = —im. The parameter s of the representation is associated to the mass scale of the system. Also 
by multiplying by — im a x (we are allowed to multiply with a x because this leaves the XXZ R matrix 
invariant) one obtains the lattice sine-Gordon L matrix [40] 

L SG (\)=( M$nK n " -im 2 sinhQuA + \ 

"" " ' ' //// 2sinli(//A ;<I>„) h ($)« ' 



h±($ n ) = l + m 2 e ±2l ^ +ili . 



LanW = 9a L S J (X) g~\ g = diag((-)^ (-)-^) (2.15) 



As was argued in [41J by taking an appropriate massless limit a lattice version of the Liouville model is 
also recovered. Indeed by multiplying the sine-Gordon L operator with 

ia ia' 

where a will be associated to the spacing of the lattice and in the classical limit a — > (see below). 
In the quantum level is a finite number and for our purposes here we set a = — i (see also section 7). 
Consider also the following limiting process [H] 

i® n — > i® n + c, fxX^fiX + c, m^O, e~ c — > oo, m 2 e~ 2c — > a 2 (2.16) 

one obtains the lattice Liouville L matrix 

(gjii„ ae — tx\— i<& n \ 

2asinh(/iA-i$ n ) h{<5> n )e- m " ) (2 ' 1?) 

h(& n ) = 1 + a 2 e- 2l ^ + ^. 

The interesting observation is that the entailed L operator (|2.17j) has a non trivial spectral (A) depen- 
dence a fact that allows the application of Bethe ansatz techniques for the derivation of the spectrum 
(see also |41j). The classical limit of the aforementioned L matrices (|2.14p . ()2. 17[) provides the corre- 
sponding classical Lax operators satisfying the zero curvature condition, and giving rise to the classical 
equations of motion of the relevant models (for more details see also |40t I4T] ) . Let us briefly review the 
connection between the quantum (lattice) versions and the classical sine-Gordon and Liouville models. 
Consider the following classical limit [4"0| I41j. the spacing a — > 0, set \i — > hfi such that ^[, ] — > {, }, and 

3 7T 3 ITT 

®n -> -x<t>{x) ~ n n ^ a-7r(x), Aim > qtti, ^lX—*u+— (2.18) 

rh is the continuum mass and 3 corresponds to the coupling constant of the sine Gordon model, and for 
the Liouville model we set ^ = 1, following the normalization of [41 j . Bearing in mind the expressions 
above (|2.18p we obtain as a — > 

L{u) = l-aU(u) +0{a 2 ) (2.19) 

then the quantities U(k) and V(u), written below provide a Lax pair for the classical counterparts of 
the lattice sine-Gordon and Liouville models. More precisely for the sine Gordon model 

i^n{x) —msmh(u + i^4>(x)) 



UW = 2Usinh(»-4* W ) 



t4>'{x) — mcosh(u + ij(j)(x)) 

rh cosh(u — i^(j)(x)) i^(j)'(x) 
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whereas for the Liouville model the Lax pair reads 

UW = d ,. 1 ' ... , V(u) = - / l / 2.21 

2 y 4sinn(« — i(p{x)) nr{x) I I \ 4cosn(u — i(p{x)j iq> (x) I 



The Lax pair satisfies the zero curvature condition 



U, V 



(2.22) 



U-V' + 

which leads to the corresponding equations of motions i.e. 

sine-Gordon model: (f>(x) — 4>"(x) H — — sm((3(f)(x)) = 



~ 2 

m 



Liouville model: <f>(x) - 4>"{x) - Aie' 2 ^ = 0. (2.23) 



Similar limiting process to (|2.16|) leads to the g-harmonic oscillator L matrix starting from 
(|2.1U|) . In fact, by simply multiplying the Liouville L matrix with an anti-diagonal matrix and bearing 
also in mind f)2. 13|) we obtain the following 

U^f""^^^ <) (2,24) 

where the operators a n , a+, V n are expressed in terms of X n , Y n as 

V n =X n , a+ = (X" 1 -qX n )Y-\ a n = Y n X n (2.25) 

and they satisfy the q harmonic oscillator algebra i.e. 

a+ a n = 1 - qVl a n a+ = 1 - q^V 2 . (2.26) 

The latter may be also seen as q deformation of the Discrete-Self- Trapping (DST) model [32]. The 
DST equation introduced in [43j to describe the non-linear dynamical behaviour of small molecules 
(ammonia, acetylene), and big molecules as well (acetanilide). Also the integrability properties of the 
DST model with two or more degrees of freedom were studied in a series of articles [341 144] . 

Finally, the simplest q deformed L matrix is obtained from (|2.14p — written in terms on the algebraic 
objects X n , Z n (see (|2.13p ) — via a simple limit 

fiX -> /xA + b, e b -> oo, - ime b -> 1 (2.27) 

then the corresponding L matrix becomes (see also [31]) 

and this L matrix was used for the construction of lattice versions of the KdV system (for more details 
on the lattice KdV see also [4311351147] ). 
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3 The open transfer matrix 



As is well known to construct the transfer matrix associated to an open spin chain one needs in addition 
to the R (L) matrix one more fundamental object, namely the K matrix, which is proportional to the 
physical boundary S matrix, acting on V, and is a solution of the reflection equation [B], 

R 12 (X 1 - A 2 ) K 1 (X 1 ) R 21 (X 1 + A 2 ) K 2 (\ 2 ) = K 2 (\ 2 ) R l2 {M + A 2 ) K^Xt) i? 2 i(Ai - A 2 ). (3.1) 



As proposed in [28], [48] an effective way of finding solutions of the reflection equations (|3.ip is to exploit 
certain algebraic structures such as (boundary) Temperley-Lieb algebras |48j-|52j. 

Let us recall how one may express the XXZ R and K matrices, solutions of the Yang-Baxter and 
reflection equations respectively, in terms of the generators of the boundary Temperley-Lieb (blob) 
algebra. First define the blob algebra bjv(q, Q), which is a quotient of the affine Hecke algebra [50J, with 
generators Ui,U 2 , ...,U/v-i and Uq, and relations: 



U,±i u z u, ±1 



5U h U Uo = 5oU 
U, ±lj UiU Ui= 7 Ui 
0, \l-k\?l 



(3.2) 



5 = —(q + q^ 1 ), q = e iAt , and Jo, 7 are constants depending on q and Q = e im ^. Physically the constants 
q and Q play the role of the bulk and boundary parameter respectively of the open spin chain, which 
will be constructed in the following. 

The generators Uj, I G {1, .., N — 1}, Uo of the blob algebra are given in the spin | XXZ represen- 
tation, i.e. let the tensor representation h : b]y(q,Q) — > End((C 2 ) 0Ar ) such that (see also [51~] ) 



h(V 



/o 




Vo 

Q q 9 
- e -Q- 1 





-Q 
1 







1 

-q~ 




\ 


/ 



(3.3) 



with h(Vi) acting non-trivially on V/ ® V/+i, (where V = C 2 ) and /i(Uo) acting on Y±. Note that 
h(Vi) I S {1, . . . , N — 1} are actually the XXZ representations of the Temperley-Lieb algebra generators 
|49| I52j. The notation here is slightly modified compared to [28]. Also, for this representation we find: 
5o = ~(Q + Q~ l ) and 7 = qQ + q^Q' 1 . 

As argued in [28] tensor representations of the blob algebra provide solutions of the reflection equa- 
tion. Hence a solution of the reflection equation (|3.1f) may be written in terms of the blob algebra 
generator /i(Uo) (|3.3p as 



(3.4) 



(h\ 2 III % III 

K\ 0) (A) = 2 sinh fj,(X - — - i() cosh fi(X- — + i() I + sinh(2^A) h(U ) 
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£ is an arbitrary constant. The XXZ spin | R matrix (|2.2p may be also written in terms of the 
representation h(\Ji), [53], i.e. 



R 12 (\) = Vu sinh //(A + i) I + sinh(/iA) h(Vi)j (3.5) 

where "P is the permutation operator, such that Via b) = 6 (g> a. 

The matrix obtained from the blob algebra (|3.4p coincides with the K matrix found in |54| [55] 
(written in the homogeneous gradation) i.e. 

R(x) = f sinhfi(-\ + iO e» x Asinlu>A) \ 

1 J y g- 9 / t sinh(2^A) sinh/i(A + ^) e^ A J 1 ' J 

subject to the following identifications 



o sinh(im/i), — — = — sinh(2i^£). (3.7) 

The -fC matrix (|3.4p . (|3.6p is written in the homogeneous gradation. To obtain the matrix in the principal 
gradation it is necessary to perform a simple gauge transformation, i.e. 

if(p)(A) = V(A) if^(A) V(A) (3.8) 

where V is given by (j2.3|) . 

Another representation that leads eventually to upper (lower) triangular K matrices follows essen- 



tially if one considers the limit q , Q —* oo, then the representation of the boundary element f|3.3|) 
reduces to: 



/>(U ) = 7 y (8) ... ® 1. (3.9) 




Indeed the matrices /i(U;) (|3.3p and /i(Uo) (|3.9p form a representation of the blob algebra with 6q = —Q 
and J = q Q. The corresponding isT matrix then is 

K(X) = ( sinh M- A + ^) e^ A g e Ksinh(2/iA) \ 
\ sinh/x(A + if) e~^ A / 

where in this case the relations among the boundary parameters become 



= -sinh(2z>C)- (3.11) 

It is clear that by considering q~ s — > oo one recovers a lower triangular matrix. In fact, such upper 
(lower) triangular K matrices may be thought of as the result of a very particular gauge transformation 
of the type (|2.15p . where the mass parameter m tx q s becomes very big (or very small). 

Finally, to construct the generating function of the conserved commuting quantities associated to 
an open spin chain, one needs two solutions X ± (A) of the reflection equation (|3.ip . Indeed, the transfer 
matrix t(X) associated to an open spin chain of N sites [7] reads 

t(X) = tr [k+(X) T (A) Kq{X) Tb(A)} = tr Q [k+(X) T (A)} (3.12) 
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where T is a solution of the reflection equation (|3.ip and 



To (A) = Lqi (A - i&)L 02 (X - *0) • • • L 7v(A - i9) To (A) = L 0JV (A + »0) • • • L 02 (A + £9)L 01 (A + i9{)3.13) 

L is one of the matrices introduced in section 2 and L(X) = L~ l {— A). As usually we keep only the index 
0, associated to the 'auxiliary' space, whereas the indices 1,...,N, corresponding to the 'quantum' 
spaces are suppressed. The parameter 9 is called inhomogeneity, and in principle one could attach a 
different one at each site of the spin chain. Also K + (X) = M K{—X — i), K is any solution of the 
reflection equation, and M is defined as: 

M = I, Principal gradation 

M = diag(q, q -1 ), Homogeneous gradation. (3-14) 

K~ is a c-number solution (|3.6[) of the reflection equation as well. The transfer matrix (I3.12p provides 
a family of commuting operators [7J, i-e. 

t(X), t(A')l = 0. (3.15) 

At this stage the quantity (|3.12|) . (|3.13|) is a purely algebraic construction. It only acquires a physical 
meaning describing a spin chain system once particular representations are chosen to act on the quantum 
spaces. 



4 Hamiltonian of a spin chain coupled to a dynamical boundary 

Already in the introduction the general form of a Hamiltonian with dynamical boundary term was 
presented (ll.lj) . We focus here in a special case of the algebraic transfer matrix (|3.12p . (|3.13j) described 
in the previous section. More precisely, we consider the situation where the spin ^ representation is 
assigned at each quantum space (L i— » R, set also 9 = 0), and the right end of the chain is coupled to 
dynamical degrees of freedom, associated to the Uq^slz), and the q harmonic oscillator algebras i.e. 

t(X) = tr {K+(X) R 01 (X) . . . R 0N (X) K (X) R 0N (X) . • • Roi(X)} (4.1) 

where R is the XXZ matrix (|2.2p . and now the K~ matrix is expressed in a factorized form [7J [26l [271 

[201 EH E3 M 

K (X) = L od (A - i0 o ) K (X) L od (A + i& ) (4.2) 

the index d characterizes the extra 'dynamical space' attached to the right boundary of the chain, K is 
a c-number solution of the reflection equation (|3.6p with parameters £~ , 9~ , L being the matrix 
given by (|2.5p or (|2.24p . and K + = MK'{—X — i) where K' is again a c-number solution of the reflection 
equation of the general type (|3.6p . but with parameters £ + , k + , + . Note that the dynamical space 
may be thought of as an extra site of the spin chain, but here it is preferable to treat it as a dynamical 
system attached to the boundary. It will be instructive for future reference to provide at this point 
explicit expressions for the dynamical matrix M appearing in the Hamiltonian (jl.ip . and due to the 
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dynamical nature of the boundary the entries of Ai in (jl.ip will be written in terms of the generators 
of the Uq{sl2) and q harmonic oscillator algebras as we will see below. In general the Hamiltonian is 

A t 



defined as proportional to -w-t(A) | a=o ; m particular we choose to normalize as (we also set = 0): 



(sinhf^))- 2 ^ 1 ^- 1 / r _i_,. A -1 d , An , 
H = - a. . ( tr o K oW) ^(A)|a=o, (4.3) 



4^sinh(i / u£~) V u 7 dX 
where T = Lq^{— iQq) Lq^zQq), and for: 

(I) the lA q {sl2) case, T = w — coshi/x(l + 20o) with w being the Casimir of U q {sl2) 

w = q A 2 + q- l D 2 + (q- q~ l ) 2 B C = q~ l A 2 + qD 2 + (q - q^fC B. (4.4) 

(II) the g-harmonic oscillator, T = 1 — g 1+2e)o . 

Then having in mind that -R(O) = i?(0) = sinhi/x V, K(0) = sinh(i^~) I and also define 

H H = -^^{Vh RuW) (4.5) 
2/j, dX 

we may rewrite the Hamiltonian as 

n_, _ ST tt sinh(i/i) T- 1 ( d ~_ m \ , tr K+(0) H 10 

H -l t Hl ^- i^mK^T) {dX K ^ X) ) A=o + tr K+(0) " ^ 

The first term in (14.61) gives rise to the bulk spin-spin interaction between first neighbours appearing 
in (II. lj) . whereas the second term gives rise to Ai, and the third one to the last three terms of (11.11) . 
We shall provide at this point the explicit values of the constants c% of (jl.ip . It is straightforward to 
compute from (|4.6p that: 



sinh(z/x) cosh(z/x£ + ) smh(i(j,)K + q e+ smh.(ifi)K + q e+ 

01 = 4sinh(^+) ' ° 2 = 4sinh(i/^+) ' C3 = 4sinh(^ + ) ' ' 

where the parameters £ + , k + , + are apparently associated to the left boundary K + . To summarize 
the Hamiltonian is finally written exactly as in (jl.ll) 



1 N ~ l N 
H = -- J2 + + cosh(i/x) of o^) - — cosh(i/i) + M f 



iV-l 

4 
i=i 

sinh(i^) cosh(i//£ + ) z sinh(i/x)K + g _6 ' + + smh{in)n + q e+ 

4sinh(^+) ^ 4sinh(i/x£+) ^ ~ 4sinh(^+) ° l ' ' ^ 

Our objective in what follows is to give explicit expressions for the matrix Ai occurring as a boundary 
contribution in the Hamiltonian above due to the second term in (14.61). 



(I) The lA q (sl2) algebra: The entries of Ai 'act' to the 'dynamical space' which in this case is asso- 
ciated to a copy of U q (sl2), hence they are expressed in terms of the U q (sl2) generators as 

sinh(i/i) J 7 ^ 1 



Mi 



4sinh(f/i£ 



cosher) fa A 2 + 2q- l D 2 -w- q 1+2e ° - q' 1 " 2 ^ 
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M 



22 



M 



12 



M 



21 



+ 2K-(q-q- l )[q- y £(g i+w ° A - q-^" D) + g y (g~ w M - q^°D)C ) + 2 sinh(i/i£~) ( ^ - q^D 2 
sinh(i/i) 



4sinh(i/x£~ 



cosh(^r) ( 2g _1 -4 2 + 2gL> 2 - w - q 1+2@0 - q~ l ~ 2 ^ 



+ 2K-(q-q- L )(q y C '(q 1+@0 D - q~ A) + g~ y (q-"°D - q"°A)B ) + 2 sinh(z/if - ) ( g£> 2 - g~M 



,-i-e 



-e . 



,e 



-1 /l 2 



sinh(i/x) J 7 



-i 



2cosh(i/^ _ )(g - g _1 )(g e °D - q~' ao A)B + 2sinh(i/x£~)(g - g'^g^'M + g w °L>).B 



-e 



4sinh(i//£~ 

+ 2/c- (g- e ~ (g - g- 1 ) 2 ^ 2 + g 9 ~ (g + (T 1 - g^D 2 - g- 1 - 20 ^ 2 ; 
sinh(i/i) 



2cosh(i/z£ _ )(9 - q~ l ){q~ e °D - q™°A)C + 2sinh(i/if~)(g " Ofa** ^ + q-"°D)C 



,e 



-e , 



4sinh(i/x£ ) 
+ 2 K -(^(g-g- 1 ) 2 C 2 +g-^(g + g- 1 



l-28o£)2 _ „l+2O ^2 



(4.9) 



(II) T/ie g harmonic oscillator algebra: In this case the entries of Ai belong to the q harmonic oscillator 
algebra (|2,26p generated by V, a, a + 



sinh(i/i) T 1 



cosh(i^ _ )(l + q l+20 ° - 2qV 2 ) + 2 sinh(z/xf-)gF 2 



22 



4sinh(i//£" 

+2«T (V e ~ +1 a+V + q e ~ (q- e °V - q^V^a 
sinh(i/i) J 7 ^ 1 



4sinh(i/i£; 



cosh(i^~)(l + g 1+2e ° - 2q~ 1 V 2 ) - 2sinh(^-)^ 1 ^ 2 



+2k 



~ 6 q e °Va+ -q°~a{q 



-l-29 



F-g 



l+20o t/-1 



v- 1 )) 



M 
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21 



sinh(i/i) 1 
4sinh(i/x£~) 

+2«~ (q 9 ~(q + q 

sinh(i/i) JF -1 
4sinh(i/x£~) 



2cosh(^-)(g _e °^ - q e °V- 1 )a + + 2 sinh(^-)(g~ e °y + g 00 ^ 1 
_ ^1+200^-2 _ g -i-2e oy2) + q ~e- {a+ f 

2cosh.{ifiC)q e °Va + 2 sinh(i^~)g 0o ya + 2«T (g e ~a 2 - g~ e ~g~ 1+20 ° V 2 ' 



(4.10) 



As expected the entries of the matrix A4, due to right dynamical boundary, are expressed in terms of 
the generators of the two main algebras i.e. Uq^sl?) and the q harmonic oscillator. These generators 
occur naturally in the expressions above due to the presence of Lq^, -^od m the dynamical boundary we 
are considering (|4.2p . To find the spectrum and the eigenstates of such a Hamiltonian one has to fix 
the representation at the boundary, and this will be done in section 6 where the cyclic representation 
(|2,12p is considered for both Uq^sl?) and the q harmonic oscillator. 



5 The boundary non-local charges 

The so called boundary non-local charges [56]-[59] will be now introduced. Such charges are going to 
play an essential role in the derivation of the spectrum as we shall see in a subsequent section. They are 
obtained via the standard procedure, by considering the asymptotics of the operator T as fiX — > 00 (see 
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also |57[|58j ). Note that in what follows the L, L and matrices are written in the principal gradation. 



(T)The XXZ and lattice sine Gordon models: The asymptotics of L(X — 0), L(X + 0) as fiX — > oo, 
bearing in mind the evaluation homomorphism (|A,5p . is given as (^O is considered finite) 

L 0n (^A - oo) oc fe \+ e -^w(^ f2 h ' n ^j+..., 

L 0n ^X ^ oc) oc ( h ' n )+e-^w( 62 > n )+... (5.1) 

with A;2 = ki . The subscript n simply denotes that the corresponding objects are associated to the n th 
site ('quantum space') of the chain, 

w = 2 sinh(^) ff~s +e ', w = 2 sinh(i/i) q^ & , ©' = + ^. (5.2) 



By also considering K , being of the form (|3.6p with boundary parameter £ , k , # , as — > oo (|3.6|) 
(the constants k~,0~ are considered to be finite) 

iT^A - oo) oc ( _ q6 K )+ e-* ( " e "^ ^ )+... (5.3) 



q k 



one may formulate the corresponding behavior of T (|3.12p . namely 

T(/iA -> oo) oc ^ ^ j + 2sinh(i/i) e""* ^ Q * ) Q (JV) ) + • • • (5-4) 

where the non-local charges Q\ N \ i £ {1, 2} are given by the following expressions |56j-|59j 

Qf = q-^K^E^ + ,He^f)^) + Xi ( K Wf, i € {1, 2}. (5.5) 

k\ N \ E^ and , provide N coproducts of the quantum Kac-Moody algebra U q (sl2) (see also 
Appendix) 

= ^ N \h), E^ = AW(e,-), = AW(/<), i € {1, 2} (5.6) 

and the constants appearing in (|5.5p are: 

x i = "o ^T7-T' x 2 = ^ r T7 - T , 0i = 0'-0-, 2 = 0' + fT. (5.7) 

2k smn(z/i) 2k smh^/xj 

It is clear that in the case of an upper triangular K~ matrix the boundary non-local charges follow 
immediately from the expressions above, given that the term proportional to q~ e is omitted. Analogous 
expressions are obtained for a lower triangular matrix, and in this case the term proportional to q® is 
omitted in the expressions (|5.5p . 



It is worth elaborating a bit further on the lattice sine-Gordon. Our aim is to provide explicit 
expressions of the lattice topological charge |40j and the fractional-spin integrals of motion for the lattice 
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sine Gordon model. To achieve this we recall that locally the L matrix (|2.14p lacks a highest (lowest) 
weight, and as a consequence the algebraic Bethe ansatz formulation cannot be applied. However, it 
was proposed in [40J that one may consider a 2N site spin chain and locally can deal with a new L 
matrix constructed as a combination of two successive L matrices. Here we are dealing with a system 
with open boundaries, and we also need to define the corresponding L matrix. Indeed we define the 
following objects 

Lon(A) = Lq 2n-l(A) Lq 2n(X), L,Qn(X) = Lq 2 n (X) Lq 2n-l(A) (5.8) 

where the subscripts 2n, 2n — 1 in the L operators denote the corresponding site in the one dimensional 
lattice. Then by keeping the lowest order terms in the /j,X — ► oo expansion the L, L matrices reduce to 
the same form as in f)5. 1 1) . In this case in particular fej, a, fi are given by the following expressions: 

U. — p -*(*2n-*2n-x) U n — 1.-1 

"-1,71 — e ) ^2,n — "'In 

ei, n = -^-^(M^nK^e^- 1 + M^n-iK^'V* 2 "), 

fl,n = ^ T (/ i+ («I>2n-i)e m2 "- 1 e i * 2 " + h_($ 2n )e- m ^e 1 *^), 

q-q 

e 2 ,n = ^ T (/ l -($2n)e- in2 "e-^- 1 + h+^-^^e^) 

q-q 1 

f2,n = ^ T (/ i+ (^2n)e iri2 "e- i * 2 "- 1 + M^n-iK^ 2 "-^* 2 "). (5.9) 

q-q 1 

In the context of the lattice sine-Gordon the N coproduct (see also Appendix A) K^ N ' plays the 
role of the lattice topological charge, whereas the coproducts e[ N \ F± play the role of the lattice 
fractional-spin integrals of motions, i.e. the lattice analogues of the classical expressions appearing in 
[60] . generating U q {sl2)- 

(II) The q harmonic oscillator. The asymptotic expansion of L(X — 0) and L(X + 0) as fiX — > oo 
gives rise to the following matrices 

L Qn (/iA -> oo) oc ^ Vn ^ +w ' e -^(^ a a "^)+..., 

L 0n ( M A oo) a ( K ^ j +w'e-^( ^ °" j + . . . (5.10) 

w ' = q-\+®\ u/ = g-f- '. (5.11) 

Finally forming the asymptotics of T, and bearing in mind (|5.10p . (|5.3p . we obtain the corresponding 
boundary non-local charges for the q harmonic oscillator 

T(/,A^oo)ocf q6 \+e~^( Qi " Q (N))+--- (5-12) 

and the explicit expressions of the boundary charges, which are new are given by, 
Ql = g-l+ei^+WyW + q -i-e lv (N) A (N) _ £^ (F (Ar) )2) Q W = q -^e 2A (N) {v (N) r i {5A3) 

AC 
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where we define 

JV-l N N 

A+W = (g) V n a+, #) = ai <g) K, F W = (g) K 

n=l n=2 n=l 

A 7 " 

^(JV) = ^i/0...0i/ 8ari8 i/- 1 ...®r 1 . (5.14) 

71=1 

Similarly if one considers an upper (lower) triangular K~ matrix then the corresponding boundary 
charges are obtained by ignoring the terms proportional to q~ 6 (q 6 ) exactly as in the previous case. 

It was shown in [58j [59] that depending on the choice of the left boundary either Q± or Q2 01 a 
combination of the two is conserved. For instance for the special case where the left boundary (principal 
gradation) is 

K+(X) = diag(e^ x+i \ e -^ A+i )) (5.15) 

the charge is a conserved quantity [58j . Actually a stronger statement was shown in [58J, i.e. 

the charge Q± turns out to be also the centralizer of the boundary Temperley-Lieb (blob) algebra 
in the spin | XXZ representation. This immediately implies the commutation of the corresponding 
Hamiltonian, written in terms of the blob algebra generators, with Q± . The derived non-local charges 
will turn out to play a crucial role in the identification of the spectrum of the corresponding models, 
hence they are not only of mathematical but of physical significance as well. Their relevance to the 
spectrum will be examined in detail in the subsequent section. 



6 Diagonalization of the transfer matrix and Bethe ansatz 
6.1 The pseudo- vacuum 

To identify the spectrum and formulate the Bethe ansatz equations of the models under consideration 
we shall apply the method employed in |19[ [33] . In the present section we deal with the XXZ model 
and the q harmonic oscillator, in the cyclic representation (|2.12p . with non-diagonal boundaries. The 
lattice sine-Gordon and Liouville models with diagonal boundaries only will be treated in a subsequent 
section. Let us give a brief account of the method, however for a more detailed analysis on the subject 
we refer the reader to [19] . One first introduces suitable gauge transformations known also as Darboux 
matrices (see e.g. 



with 



M n (A) = (x n (X), y n (A)), M„(A) = (x n+1 (X), 5^i(A)) (6.1) 

*n(A) = y 6 ^ J , Y n (X) = ( 6 ^ Vn V x n = x e-^ n , y n = y e^ n . (6.2) 

Here {xo,yo} £ C are A— independent and will be determined explicitly in the following. Introduce also 
the matrices 

( ) • ( t:« ) ■ <6 - 3) 
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with 



X n (X) = — ( - e»\ x n ) , Y n (X) = —L- (>\ - y, 




■ n ] i 

Xn Un ^ ' X n Vn 

X n (\) = ( - e"\ x n ) , Y n (X) = (e"\ - y n ) . (6.4) 

X n - Vn-2 \ ' X n+ 2 ~ Un ^ ' 

One can check that certain combinations of (|6.2p . (|6.4p lead to orthogonal relations among the various 
vectors. Also, face- vertex correspondence relations corresponding to commutation relations between the 
above quantities can be derived using the R— matrix (|2.2p . which are omitted here for brevity [19J. We 
now introduce the gauge transformed L— operator, 

L n (m\X) = M-^_ l)g+m {\)L n {\)M ng+m {\) 

S n (m\X) = M- 1 +m (-A)L- 1 (-A)M (n _ 1)s+m (-A)= ( f) . (6.5) 

where g depends on the choice of representation. A priori one may associate a different representation 
at each site of the chain, here however we consider for the sake of simplicity the case where all sites 
correspond to the same representation. The occurrence of a different g at each site, depending on the 
choice of representation, is a non trivial observation facilitating the generalization of the formulation 
described [19], where all sites are associated to the fundamental representation. Note for instance, that 
for the spin | representation the factor g is unit, for the spin s representation is 2s, and for the models we 
are considering here in the cyclic representation will be specified later on. Also at the end of the section 
6.2 we shall discuss the special case where all sites correspond to the spin 5 and the space attached to 
the boundary is associated to the cyclic representation (|2.12j) . This scenario in fact corresponds to the 
Hamiltonian described in section 4 



t{\)=tr [k+{\)%(\)Y (6.6) 



It follows that the transfer matrix (|3,12p can be rewritten with the help of the aforementioned gauge 
transformations (|6.5p . and related orthogonal relations [19] as 



where 
K + (m\X) 

and 




Y m (-\)K+(\)X m (\) Y m (-\)K+{\)Y m -2{\) 
X m (-\)K+(\)X m+2 (\) X m (-\)K+(\)Y m (\) 



(6.7) 



f(X)= AmW BmW I = I ^m- 2 (A)T(A)X m (-A) Y m (\)T(\)Y m {-\) \ 
\ C m {\) V m {\) ) \ X m {\)T{\)X m {-\) X m+2 (\)T{\)Y m (-\) ] ' 

Similarly to [19j . one defines the transformed K~ matrix as 

k-(l\X)= ( \ ( Yl-2WK-(X)X t (-X) Y l (X)K-(X)Y l (-X) \ 

1 1 ' { K^(l\X) K 4 {l\X) j I X l (X)K-(X)X l (-X) X l+2 (X)K-(X)Y l (-X) I 1 ' 1 
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with I = m + Ng. In the case where a different representation is assigned at each site then I = 
171 + Y2 n =i9m wnere 9n characterizes the representation at each site. Note that here the matrices 
are c- number solutions of the reflection equation, with relevant parameters k^, q 9± . 

The first step for the diagonalization of the transfer matrix (|3.12|) is the construction of an appro- 
priate pseudo- vacuum. Such a state will be of the general form: 

N 

nH = w w (6.10) 

n=l 

Here we denote Un the local pseudo-vacuum annihilated by the lower left elements of the transformed 



matrices (j6.5j) . i.e. 

%, f n 4 m) =0. (6.11) 

Such a choice ensures that the operators T(A), T(A) (I3.13p . acting on the pseudo- vacuum state reduce to 
upper triangular matrices. The above requirements (|6.1ip lead to certain sets of algebraic constraints. 
In particular, let 

L 0n (X) = ( MX) Pn{X) ) (6.12) 
{ ] \ 7n(A) 8 n (\) J { ' 

from the action of 7 n on the local pseudo-vacuum we obtain: 

- Xm+ng+i a n {\) + x m+(n _i )g+1 S n (X) + e~ fMX x m+ng+1 x m+(n _i )g+1 7„(A) - e^ x /3 n (X) u^ m) = Qj6.13) 

from the action of 7^ a similar constraint is entailed. 

The non-diagonal elements of (|6.7|) . (|6.9|) acting on the pseudo- vacuum state (|6.1U|) must satisfy 



certain constraints such that the pseudo-vacuum is indeed an eigenstate of the transfer matrix 
Indeed, we associate the integers m° and ijiq to the left and right boundaries respectively, and impose 
the following: 

K+(m°\X) = K+(m°\X) = , k^(mo\X) = , (6.14) 

K^ 3 are given in (|6.7|) . (|6.9|) . The above constraints may be solved [TO] fixing the relations among the 
parameters mo, m°, xo,yo- Indeed let xq = — ie~ ifJ, ^ +rr >, yo = — ie _i/ ^~ 7 ), then the relations among 
the boundary parameters (|6.14p are 

-i coshz^(/3 + + 7 + + 1), = —% coshi^(/3 + — 7 + + 1), 



2k+ ^ 1 n 2rH 

i cosh ip,{p~ + 7~), (6.15) 



2k 



and we define 



For relevant discussions on these constraints see also |20| , I18 |, [2T | ll9j. It would be illuminating to consider 
the above constraints bearing in mind the identifications of the boundary parameter with the parameters 
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of the blob algebra [58]. Indeed analogous relations emerge in the double Temperley-Lieb algebra for 
exceptional values of the algebra parameters in [23] . It is not clear to us however why such constraints 
hold true for generic values of the boundary parameters. The natural question raised is whether this 
is simply a disadvantage of the methods employed or it reflects a deeper physical or algebraic reason. 
This is an intricate point, which merits further investigation. 

In what follows we examine the pseudo- vacua for both the XXZ model, and the q harmonic oscillator. 
Their explicit forms are specified by solving the formulas (|6.13p . which lead to sets of difference equations. 

(I) The XXZ model: As already mentioned we shall treat the cyclic representation (|2.5p , (|2.10p , (|2.12l) , 
for the special case where q is root of unity. We associate each site of the chain with the cyclic repre- 
sentation, and we define the local pseudo-vacuum state as 



1=1 



/, is a p dimensional column vector with zeroes everywhere apart from the I th position. The constraints 
(|6.13p provide (a) the value of g = 2s — 2 and (b) recursion relations among the u>£ m : 

w fc+'i n) _ g s ~ 1+e sinh jjjtfs - s + 1) wj" 1 '^ _ -o+i~ s sinh(i/is) 

W M Xm+[n _ 1)g+l S i n h^(fc + S ) ' W M x m+(n-l) 9 +l 1 ginh _ s) V>-**) 



By normalizing the state with w^" 1 '^ = 1 we can write 



w k 



(m,n) _ ( (f 1+0 \ k ~ x TT smhi/x(j -5 + 1 



Xm+{ n -i)g+i 1 fj[ sinhi^(j + s 



n r", v • (6.i9) 

A -L (jinn 9/(11 -I- el 



Comparing the expressions of w^'^ from equations (|6.18p , (|6.19p we obtain the following constraint on 



( 



x m+(n _ 1)g+1 q- e+1 - s y = 1. (6.20) 
Relations (|6.18|) are sufficient for the derivation of the spectrum. 

(II) The q harmonic oscillator. For the L matrix describing the q harmonic oscillator (|2.24|) we consider 
X — > q z X, Y — > q%Y (cyclic representation ()2.12p ). In this case g = —1. The local pseudo-vacuum will 
be also of the form (|6.17p . and the corresponding recursion relations are also provided by solving the 
constraints (|6.13j) . 



w 



2sinh^(fc-2 + ~) w^' n) __ x m+ng+1 q z 1 H 



fc+i 



w ^ n) ~ V»s+i f Hz " 1_e ' 4 m ' n) " 2sinh^(i-z) 
Again normalizing the pseudo- vacuum with u^ m ' n ' ) = 1 we conclude that 



(6.21) 



' ;i) - ' ) sinh //,(./ z + \), (6.22) 
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and the constraint necessary from cyclicity leads to: 



1 ^ 1 

-x m+1+ng q z - l ~ e ) P = JJ q k sinh ifx(k+--z). (6.23) 



2 — u <,~r ±~r i i>y i. j XX ' ^ 

k=i 

To derive the spectrum of the transfer matrix for both the XXZ and the q harmonic oscillator we 
also need to know the action of the transformed diagonal elements on the pseudo-vacuum, which are 
given by: 

4"0 4™) = a(A) 4 m+1 \ 

tfM W M = X m+ng+l ~ Vm+ng-1 ^ 

•^m+(n—l)g+l Um+(n— l)g— 1 

c# TO > a'(A) 4— V, 

~ 5 (>m) Jm) = XmHn-l)9-ym+(n-l)g ^ ^m+l)_ ( g ^ 

S-m+ng Um+ng 

The values of a, b, a' , b' for each representation are given below 

(I) a(A) = q~ 3+1 sinh //(A + is - iQ - i), 6(A) = q 3 " 1 sinh fi(X - is - iQ + i) 
a' (A) = gT 1+s sinh/i(A + is + iQ), b'(X) = q 1 ' 3 sinh / u(A - is + iQ + 2i) 
(II) a(A) = igi-^+Qe-M, 6(A) = g 2 ~^ sinh/i(A - iQ + ^) 

a'(\) = q-^ +z smhn(\ + ^ + iQ), b'{\) = p 2 ~ z e" {x+ie) . (6.25) 
6.2 Spectrum and Bethe ansatz equations 

Now that the appropriate gauge transformations and the corresponding pseudo-vacua have been explic- 
itly derived, we can proceed with the computation of the eigenvalues of the transfer matrix. Namely, 
we are looking for the solution of the following eigenvalue problem 

t(A) * = A(A) ^ (6.26) 

where \E' is the general Bethe ansatz state of the form 

* = Bm°-2(Ai) • • • S m o_ 2M (A M ) ^ (m) , (6.27) 

and m will be defined below [19]. Then one has to solve a typical eigenvalue problem (|6.26p written 
explicitly as 

t(X) f = (ir+(m°|A)Ano(A) + ^ 4 + (m°|A)P m o(A)) *, (6.28) 

where 

~ sinh u(2\ + i) sinh u(im + i^f + i) „ sinh u(im + i7 + 2A + i) , 

2?m(A = , r . r , , ' L V m (X) \ - , ' \ ~An(A) 

smn^j sinni J u(m + 7j smJii/i(m + 7) 

and the elements ifj 4 (m°|A) are given in Appendix B. Using the definition of T(\) in (|3. 12|) and ()6.8j) . 
the action of the operators A m (X) and T> m (\) on the bulk part of the pseudo- vacuum (|6,10p may be 
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derived. Assuming that each site n of the chain is associated to either the XXZ model (|2.5p or to the q 
harmonic oscillator (|2.24p in the cyclic representation (|2,12p . it is convenient to define 

f n (A) = a(A) a'(A), h n (A) = b(X) b'(X), (6.29) 

a, b, a', b' are defined in (|6.25p . Then bearing in mind relations (|6.24p and appropriate algebraic 
relations arising from the reflection equation (see also [19]) we obtain: 

N 

A m {\)^ m ) = l[f n (\) K-(m \\)rt^ 
n=l 
N 

V m {X)^ = l[K(X) K7(m \X) rt m \ (6.30) 

n=l 

mo = m + gN, and iT-^ (mo | A) are introduced in Appendix A. 

Having in mind the structure of the general Bethe state (I6.27P it is obvious that we need exchange 
relations between the operators A m (X) and V m (X) with £>;(A), which can be deduced with the help of 
the reflection equation (|3.ip (see also [E]), then one obtains 

, , sinh^(Ai + A 2 ) sinh^(Ai - A 2 - i) f . , . 
A m+2 {X l )B m {X 2 ) - sinh ^ (Ai _ Aa) sinMAi + A2 + i) B m (A 2 )A»(Ai) 

sinh(2/iA 2 ) sinh(i^) sinhyu(Ai - A 2 - — 17 - i) , . , . 

£> m (Al JAn (A 2 J 



sinh/i(zm + Z7 + i) sinh/i(Ai — A 2 ) sinh/i(2A 2 + i) 
sinh(i/i) sinh(i/i) sinh/i(— Ai — A 2 + im + 27) 
smh. fi(im + ij + i) sinh/i(Ai + A 2 + i) sinh^(2A 2 + i) 



^ m (Ai)P m (A 2 ) (6.31) 



and 



-n t\ \k t\ \ sinh/i(Ai + A 2 + 2z) sinh/i(Ai - A 2 + t) ~ 
P m+ 2(A!)^ m (A 2 ) = sinh/l(Ai _ Aa) sinh ^ (Ai + A2 + i} H m (A 2 )^ m (A 1 ) 

sinh//(2Ai + 2i) sinh(i/x) sinh / u(Ai - A 2 + im + £7 + i)^ ^ 

sinh /i(im + 17 + i) sinh^(Ai — A 2 ) sinh^(2A 2 + i) m m 

sinh(2//A 2 ) sinh 2//(Ai + i) sinhyu(Ai + A 2 + im + i-y + 2%) . . . . 

H ■ u r T~ — I -\ T~ / \ 1 \ — , -x r-r , . — -tt— O m (Ai)^ m (A 2 ) . (6.32) 

smn + 17 + 1) smn/i(Ai + A 2 + z) smn^(2A 2 +z) 

In fact all the above exchange relations, arising from the reflection equation, hold exactly as described 
in [19]. What is only modified is the action of the diagonal elements and consequently of A m , T^m on 
the pseudo- vacuum (|6.30p . which clearly depends on the choice of the particular representation. It may 
be shown that the state ^ is indeed an eigenstate of the transfer matrix if we impose m = m° — 2M. 
Without loosing generality we assume for simplicity m° = 0, and then it immediately follows that 

m = Ng -2M . (6.33) 

We may now deduce the transfer matrix eigenvalues by virtue of equations (|6.28p . (|6.30p . (|6.3ip . (|6.32p 
i.e. 

. / \ \ / T ^/ , . . t~t /,x tt sinh u(X + A,) sinh a(X — A,- — i) 

A A = Kt (O A) K7 (m A) IT f n A) 1 T — -?-)- — — — - 3 

v ' \ 1 v 1 ' 1 v 1 ; y y 1 11 sinh/i(A + Xj + i) sinh / u(A-A j ) 

71 — 1 j — 1 

T , . . , . _ , , , . T-r , /A . J^- sinh u(A + A 7 - + 2i) sinh u(A — A 7 - + i) \ 
+ At(0[A)AV(m |A)nK(A) g a^jpj?^ sinh ,(A-A,) ) ( ° 4) 
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provided a certain combination of 'unwanted' terms arising from the commutation relations (|6.3ip . 
(|6.32p that appear in the eigenvalue expression, is vanishing. This is true as long as Aj's satisfy a set 
of conditions, namely the Bethe ansatz equations. These equations guarantee also analyticity of the 
eigenvalues, and they are written in the familiar form 

N f (\ — il M 

H(Xi) n i n _\\ = - n - w + ^ > ( e - 35 ) 

n=1 ^n{\ 2 > j=l 

where we define 

V ; sinhyu(A-f) V ' K 4 + (0|A- t)K-(m \\- §) v ' 

A few comments are in order at this point. Notice that for the q harmonic oscillator a factor q~ N , which 
reduces to ±1 for N = kp (k integer or half integer), occurs in the left hand side of the Bethe ansatz 
equations, resembling the case of the XXZ chain with twisted boundary conditions. As realized in [UJ 
a similar factor appears also in the Bethe ansatz equations of the lattice Liouville model with periodic 
boundary conditions. Furthermore, for particular values of the parameter s and the inhomogeneities 
i.e. (I) s = 1, (II) G = (see also (|6.25p ). the Bethe ansatz equations become degenerate, and as a 
consequence the Bethe ansatz formulation is not appropriate any more for deriving the spectrum. In 
this case one has to deal with the relevant Baxter operators and the Baxter equations [I] in order to 
examine the associated spectra. 

The case of dynamical boundaries studied in a series of papers (see e.g. [3 [T3], and 
references therein) can be regarded as a special case of the above description. Indeed, let us consider 
a typical example of dynamical boundaries, which is described by the Hamiltonian of (|4.8p . Consider 
a chain with N + 1 sites, and assign to all quantum spaces n G {1, . . . , N} the spin ^ representation, 
whereas the N + 1 site is associated to the cyclic representation, of U q (sl2) or the q harmonic oscillator 
algebra that is we consider a dynamical boundary at the right end (see also section 4). In fact, the 
N + 1 space is the 'dynamical' space d of section 4. Then the spectrum takes the form (|6.34p . but with 
the functions f n , h n being (set here for simplicity all the inhomogeneities to zero): 

f n (A-^) = sinh 2 / /(A + ^), h n (A-^) = sinh 2 /i(A-^), ne{l,...,N} 

WA - ~) = a(A - I) a'(A - I), 1WA - \) = 6(A - \) b'(X - I), (6.37) 

and recall a(A), a'(A), 6(A), b'(X) are given in (|6.25p for both U q (sl2) and the q harmonic oscillator. 
The pseudovacuum for the spin | case have been found in [19] so we do not repeat them here. The local 
pseodovacuum for the site n is given by (|6,17p and (|6.19p or (|6.22p depending on the dynamics at the 
boundary i.e. Uq^sl?) or q harmonic oscillator respectively. The general Bethe states are given in (|6.27p . 
We could have considered a dynamical boundary in the left end as well, then fi = fjv+i, hi = hjv+i- 

6.3 Asymptotics and derivation of M 

In |18(. [T9] the spin ^ XXZ model with both boundaries being non-diagonal is discussed in detail. In this 
case via the asymptotic behavior of the transfer matrix it was shown that the value of the integer M is 
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fixed (see also |20l I21j for the spin s representation). Let us now consider a special left boundary given 
by (|5,15p . then the value of M may be derived in terms of the spectrum of the conserved quantity Q±*\ 
Indeed the asymptotic expansion of the transfer matrix (|6.34[) as /j,X — > oo, taking also into account 
(|B.6|) . is given for (I) the XXZ model, and (II) the q harmonic oscillator: 

(I) A(/xA -»■ OO) = -i K~ e ^KN+l)+iKN+l) ^ e iti(f3-+-y--2(s-l)N+2M) + e -irfp-+ J --2(a-l)N+2M)^ 

(II) A(/xA -» OO) = -i K~ e 2M(iV+l)+i M (27V+l) e i Ai (/3-+ 7 -+2M) (g gg) 

comparing with the results of section 5, and taking into account appropriate normalizations for the K 
and R matrices, we conclude, 

(I) Q[ N) = i Shlh- 1 ^) (y-+7-+2Af-2( S -l)iV + q -p-- 7 —2M+2(s-l)N^ 
(II) Qf ] = 2 i /-+7-+2M+iV_ ( g 39) 

The expressions of the boundary charges for both XXZ and q harmonic oscillator are given by the 
expressions (|5.5p and (|5.14p respectively. Also the sum j3~ + r y~ is related to the right boundary 
parameters of the chain (|6.15p . From equation (|6.39p the upper (lower) bounds of M (integer) are 
identified by means of the spectrum of the non-local operator which plays the role of q s * , the 

corresponding conserved quantity when diagonal boundaries are applied in the spin chain. The relation 
between and M underlines the significance of the algebraic object Q± into the derivation of the 
spectrum of the open spin chain with special boundary conditions. The identification of the spectrum 
of the boundary non-local charge is an intriguing problem, and some progress has been already achieved 
towards this direction, for particular representations in [23], however the problem in its full generality 
remains open. In fact, the operator Qi locally consists of tridiagonal form matrices for the spin s 
representation (|2.7p (we provide an example in Appendix C). The explicit diagonalization of Q± for 
generic number of sites will be discussed elsewhere. 

Similar algebraic objects, forming the corresponding boundary quantum algebra [ST], are necessary 
for the identification of the spectrum of models associated to higher rank algebras with non diagonal 
right boundary. 



7 Boundary Lattice sine— Gordon and Liouville models 

In this section we give a flavor on the lattice sine-Gordon and Liouville theories with diagonal boundaries, 
that is in the expressions for the c-number matrices we set = 0. In particular, we identify the 
spectrum and the corresponding Bethe ansatz equations for both models. 

In what follows we consider for simplicity the inhomogeneity G to be zero. We find local pseudo- 
vacua annihilating the lower left element of L, L, being effectively 'highest weight' states. Let such a 
state be of the form [lOj 0T] 

W« = /(*2n) <5($2n - $2n-l ~ n) (7.1) 
7n, In W„ = (7.2) 
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and consequently /($) should be a solution of the following set of difference equations 

h-($) 

Lattice sine-Gordon: /($ + y) = — - 7^y /W 

Lattice Liouville: /($ + /i) = (7.3) 

which occur exactly as in the bulk case. When q is root of unity a cyclicity (periodicity) requirement is 
necessary to be imposed on the pseudo- vacuum as well, i.e. (see also [10JHT]), 

/(*+PM) = /(*)• (7-4) 

To find the spectrum of the lattice systems we also need the explicit action of the diagonal entries of 
the L, L matrices on the local pseudo- vacuum. Indeed it is straightforward to show that 



where 





(A) uj n = 


= a(X) to n , ot n (X) to n = a (A) u n , 




S n (X) uj n = 


^ 6(A) w n , 5 n (A) w n = 6'(A) u n 




Lattice sine-Gordon: 


o(A) = 


= —4m 2 cosh //(A — — + — ) cosh /i(A — — — 


ir 
~2* 




6(A) = 


- —4m 2 cosh //(A + - + — ) cosh y(X + — — 


ir 
~2' 




a'(X) 


% IT* % 

= —4m 2 cosh y(X + — + — ) cosh //(A + - - 


ir 




b'(X) 


= — 4m 2 cosh/i(A + — + — ) cosh/i(A + — 


ir 


Lattice Liouville: 


a(A): 


= _ e -M(A-|) sinh/U (A- 1) 






6(A): 


= _ e -MA+f) s i n h^(A + |) 






a'(A) 


= e MA+f) sin h^(A + ^) 






6'(A) 


= e M(A+f) sinhAi (A+|), 





(7.5) 



(7.6) 

the constant r appearing in the sine-Gordon eigenvalues is given by 2coshi/ir = (m 2 + m -2 ), it plays 
the role of inhomogeneity and it is associated to the mass scale of the system. The standard algebraic 
Bethe ansatz may be applied [7j, and the spectrum is obtained having the familiar form of (I6.34h 
(K^ 4 (mo\X) — > i£~j~ 4 (£~|A), K^~ 4 (0\X) — > ifj 4 (^ + |A)), while the corresponding Bethe ansazt equations 
as customary follow as analyticity conditions on the spectrum and they have the form (|6.35j) . (16, 36ft 
with 

f n (A) = a(A) a'(A), h n (A) = 6(A) 6'(A) 

^-(r|A) |A) = ^g±^sinhM-A + zr) sinh^A + z^) 

K 4 -(r|A) A" 4 + (C + |A) = si ^ ( 2A + i) Sinh/i(A + ir + * } Sinh '"(- A + ^ + " 0. ( 7 - 7 ) 
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and a(A), a'(A), 6(A), 6'(A) given by (|7.6p for both the lattice sine-Gordon and Liouville models. The 
integer M appearing in the Bethe ansatz (|6.35p this time is associated to the spectrum of the topological 
charge k[ N) (UlSD, (JHHD for the sine-Gordon model. It should be stressed that the identification of the 
boundary non-local charges for the case of the open lattice Liouville theory is an intriguing task, which 
needs to be further pursued. 

In the same spirit the case where the right boundary is upper (lower) triangular may be examined. 
An obvious reference state still exists, and the Bethe ansatz equations may be easily extracted being of 
the familiar form. In this case of course the corresponding conserved quantity for the lattice sine-Gordon 
will be the charge Qi (|5.5p (for a trivial left boundary (|5.15|) ). but with the terms proportional to 
Q {(f ) omitted. Finally the more general scenario with non-diagonal boundaries may be studied 
employing appropriate local gauge transformations. The spectrum is anticipated to be of the form (|6.34l) . 
but this time the integer M should be associated to the corresponding conserved quantity. In general 
M being associated to the corresponding conserved non-local charge reflects essentially the underlying 
symmetry. Moreover, in this case the derivation of the reference state involves the solution of more 
complicated sets of difference equations, and it will be left for future investigations. 

8 Conclusions 

The main aim of the present study was the investigation of a class of open quantum integrable models 
emerging from the generalized XXZ model (12. 5p . In this spirit a particular example of a spin ^ XXZ 
Hamiltonian coupled to a dynamical boundary (|4.8I) was first considered. Explicit expressions of the 
corresponding dynamical boundary term associated to the two main algebras i.e the U q {sl2) and the 
q harmonic oscillator algebra were derived (|4.9p . (|4.10p . The main point is that the expressions of 
the boundary terms (14.9P and (14,101) are generic at this stage that is independent of the choice of 
representation, hence they can be exploited for a variety of dynamical systems (quantum impurities) 
coupled at the end if the chain. As a step towards the identification of the spectrum of the various 
models under consideration, and after recalling known expressions of the boundary non-local charges 
for the generalized XXZ model, we derived novel expressions of boundary non-local charges for the 
lattice sine Gordon model (|5,5p and (|5,9p . and for the q harmonic oscillator (|5.13p and (|5.14p . The 
significance of such quantities and their relevance to the spectrum is emphasized in section 6.3. 

In order to examine the spectra and eigenstates of the XXZ model and the q harmonic oscillator 
with non diagonal boundaries in the cyclic representation (|2.12p . we generalized the approach presented 
in [19]. The crucial observation is that although the local gauge transformations (see (|6,5p ) act on the 
'auxiliary' space — being always associated to the fundamental representation oiU q {sl2) — they depend 
explicitly on the choice of representation for each quantum space (site). In particular, there exists a 
parameter g incorporated in the local transformations (|6.5p . whose value depends clearly on the choice 
of representation for the each site (for the spin ^ case is unit). Based on this observation we were able 
to derive explicit expressions for the pseudovacua (|6.17j) . (|6.19p . (|6.20p . (|6.22p . (|6.23p by solving sets 
of involved recursion relations (|6.18p . (|6.2ip . In fact both the value of g and the recursion relations 
(|6.18p . (|6.2ip were deduced from the solution of the constraint (|6.13p . which is the fundamental object 
in this construction. More precisely, the main requirement imposed was that the transformed L matrices 
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(|6.5p become upper triangular after acting on the suitable local pseudo vacuum, which eventually led 
to the constraint f)6. 13|) . Consequently the Bethe states (|6.27p were identified, and the spectra (|6.25|) . 
(|6.34p . (|6.29p and Bethe ansatz equations (|6.35p of the aforementioned models were derived. The 
derivation of the pseodovacua and Bethe states is of particular interest especially in the context of 
computing correlations functions (see e.g. [H]). Moreover in the semiclassical limit as known Gaudin 
type Hamiltonians arise, whose eigenstates are the semiclassical limits of the Bethe states, and they 
satisfy Knizhnik-Zamolodchikov equations (see e.g [62]). It is a compelling task to identify the exact 
type of the Knizhnik-Zamolodchikov equations for the semiclassical limit of the models investigated here, 
this however will be undertaken elsewhere. We should also point out that within the spirit described 
above the Hamiltonian (|4.8I) . with a fixed representation (cyclic) at the right end, was treated as a 
special case and the corresponding spectrum was presented in (|6.34p . (|6.37p . 

Another intriguing new result is the link between the spectrum of the open transfer matrix with a 
trivial left boundary and a general non diagonal right boundary, and the conserved quantity (see 
(|6.38p . (|6.39p ). which plays essentially the role that S z plays in models with both diagonal boundaries. 
The study of its spectrum was also emphasized and examples on the diagonalization of the one site charge 
for the spin s and the cyclic representation (|2.12p of the U q (sl2) were also presented (see Appendix C). 
Finally, we were able to identify the spectrum (j6.34j) and (|7.7p and the corresponding Bethe ansatz 
equations (I6.35|) . (|7.7p for the lattice versions of the sine Gordon and Liouville models, but only with 
diagonal boundaries. We also realized that the corresponding pseudovacua (|7.ip . (I7.3p . (|7.4p have the 
same structure as in the periodic case |401 141] . Note that this is the first time that these models 
with open boundaries are examined, so the aforementioned results are quite useful especially from the 
point of view of computing the corresponding exact boundary S matrices and boundary thermodynamic 
properties. 

It is worth emphasizing that all the models under consideration share a common spectrum form 
(I6.34p . This is somehow anticipated given that they all belong to the same universality class, emerging 
from the generalized XXZ model (12. 5jl . More precisely, the terms JL=i m A6,34p . and consequently the 
right hand side of the Bethe asnatz equations (|6.35p . are generic that is independent of the choice of 
representation. On the other hand the terms of the spectrum including f n ,h n and the left hand side 
of the Bethe ansatz equations depend clearly on the choice of representation, hence the various models 
naturally give rise to distinct expressions. Indeed notice the various f n ,h n functions for each model 
(|HT25|> . (|Q9]h (jOTp . CLU), dZTT]). Furthermore, the spectrum AO!) , and the left hand side of the Bethe 
ansatz equations (|6.35p . depend on the choice of boundaries see e.g. in (I6.34p . and also the function 
H(X) (|6,36p . (|7.7p associated to the corresponding boundary conditions. It should be finally stressed 
that in order to obtain 'numbers' for the spectrum one has to solve the corresponding Bethe ansatz 
equations. This can be achieved numerically for finite size chains, or using thermodynamic techniques 
as N — > oo, obtaining consequential (boundary) scattering information and (boundary) thermodynamic 
quantities (see e.g. [20]). This is a particularly interesting aspect, but it is beyond the intended scope 
of the present study, and it will be pursued in a forthcoming work. 
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illuminating discussions. I wish also to thank the organizers of the 3rd Annual Network Meeting, 'Inte- 
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A The quantum Kac— Moody algebra l( q (sl<2) 

We briefly review some basic definitions concerning the quantum group structures. Let 

<■*)-( _ 2 2 - 2 2 ) (A.1) 

be the Cartan matrix of the affine Lie algebra sfa [63] , and also define 

= ^2 (A.2) 
q-q 1 

Recall that the quantum affine enveloping algebra U q (sl2) = A has the Chevalley-Serre generators 
[531 El] e «> fit h, i £ {1, 2} obeying the defining relations 

hi kj = kj ki, ki &j = q% *^bj k^, k{ fj = q ^ lJ fj ki, 

~ei, fj] = Sij k *~^* , 2} (A.3) 

and the q deformed Serre relations 

Xi Xj ~ [3] 9 Xi Xj Xi + [3] g Xi Xj Xi ~ Xj Xi =0, Xi S {e», /»}, i / j. (A.4) 
There exists a homomorphism called the evaluation homomorphism |53j n\ : lA q (sl2) — ► M q (sl2) 

7TA(ei) = ei, 7Ta(/i) = /i, 7TA(fcl) = fcl 

7r A (e 2 ) = e-^c/i, vr A (/ 2 ) = e^W, 7r A (fc 2 ) = fcf 1 , (A.5) 

c is a constant. As mentioned this algebra is also equipped with a coproduct A : .4 — > .A <8 .4, in 
particular the generators form the following coproducts 

A(y) = h®y + y®kr l , y e {*, fc} and A(kf 1 ) = kf l ®kf l . (A.6) 

The L-fold coproduct is derived using the recursion relations 

A< L > = (id® A^-^A (A.7) 

and thus one may obtain explicit expressions of the L coproduct as 

L L 

& {L) {y) = ^k l ®...®k i ® y sfcr 1 $ . . . ® y e { ei , /<}, A^A*) = (A.8) 

n=l ... n=l 

n position 
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B Transformed K matrices 



Explicit expressions for diagonal elements of the transformed K matrices are provided below 



K7(m u \X) = — — — sinhiufl -7" 1 ") H — sinh fi(2X + ij + + i) + 2i cosh(z/j/5 + ) sinh(2^A) 

sinn(i^7 + ) V 2k + 2k + 

A^"(m°|A) = — — — -^1— — sinhi/x(l +7+) + sinh/i(2A - 27+ + i) + 2i cosh(i///3 + ) sinh(2/iA) 



*r(mo|A) = sinhi/i(7 - _ i) + 7-) + sin h/ x(2A - i 7 " + i) 

— 2i cosh + 1) sinh(2^A)^ 

K 4 -(m |A) = ginh ^p+i) Up Sinh + 7 " } " IF Slnh M2A + * 7 ~ + l) 

+ 2t cosh ifi(p- + 1) sinh(2/xA) V (B.2) 



By restricting our attention to the case where the left boundary is a trivial diagonal matrix (|5.15p , then 
the transformed diagonal entries become 

Kt(m°\\) = — simWl - 7 +), ET+(m°|A) = ^ -t sinhz/x(l + 7+). (B.3) 

1V 17 2sinh(^7+) PV ' ; ' 4V 1 ' 2sinh(i/x7+) PV ' ; v ; 

For the derivation of the spectrum we shall also need certain combinations of the diagonal transformed 
elements i.e., 

ir+r \\\ ts+( \\\ , sinh//(i 7 + m + 2A + i)sinh(z/x) ~ + . 

a, m A = K7(m\X) H : ; r : — 7— -Kj(m A) , 

1V 1 ' sinh + 772 + 7) sinh fi(2X + i) 4V 1 ; ' 

K +( m \ X ) - sinh ijxjjn + 7) sinh(i/i) - + 

^ 4(m|A) " sinh^(l + m + 7) smh^(2A + i) K * (m|A) ' (R4) 



iff (mo I A) = Ki(m \\) 
sinh ifi(mc 
sinhi 

sinh /J,(imo + ij + 2A + i 



sinh z/j(m + 7 + 1) sinh /i(2A + t) 7 _ 

K 4 (mo A) = . . . JT 4 (mo A) 

8innj^(mo + 7) sinn^j 



AT (m 1 A) . (B.5) 



sinhi/x(mo + 7) 

It is instructive for our purposes here to consider the asymptotic behavior of K^. We shall only deal 
with the case where K + is a trivial diagonal matrix (|5.15p . Considering the case where fj,X — > 00 we 
conclude that: 

Kr(mo\n\ -» 00) = -i K - e 3 ^ A -W-+7-) ) K7(m |MA -» 00) = ^_^ e 5M+^(/3-+ 7 -+2) 

2 smh(i/ij 

if^mVA -» 00) = e ~^ x+i ^, K+(m°\n\ -» 00) = 2e- 3 ^ A ~^ sinh(^). (B.6) 
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C The boundary charge as a tridiagonal matrix: diagonalization 



In this appendix we consider a particular example of the boundary operator Q\ (|5.5p . We deal with its 
spin s representation (|2.7p and for only one site. In this case the operator reduces to an n x n tridiagonal 
matrix (see also |35j ) whose diagonalization will be the main objective in this Appendix 



(A) Let us first present the boundary operator for one site Q\ (|5.5p in the spin s representation ([2,71) , 
Let us set for simplicity 



A fc = q 2 



Ck, Bk 



the boundary operator may be written as 

n-l 



Qi = ^2 Ak efc fc+i + ^ B k Ck k-i + ^2 ^ k Bkk 



(C.l) 



(C.2) 



k=\ 



k=2 



k=l 



which is indeed a tridiagonal matrix. Let us now solve the corresponding eigenvalue problem. Let 
\P = Ysi=i w i fi be an eigenstate of Q\ then 



Qi * = e * 



(C.3) 



where e is the corresponding eigenvalue. The latter equations can be written in a more explicit form as 
a tridiagonal (Jacobi) matrix, 

/ Ci A\ \ 

B2 C2 Ai 



Bk Ck A k 



\ 



B n -i C n _i A n -x 

B n C n J 



( W\ \ 

W 2 



\ w n / 



( w l \ 

W2 



(C.4) 



Note that in [35j the diagonalization of triadiagonal matrices and their relation to q hypergeometric 
series is discussed (also related to Leonard pairs, Askey- Wilson polynomials etc. see e.g. [59] and 
references therein). 

One has now to solve the following recursion relations, in order to determine the factors W[ of the 
eigenstate, 

Ai u> 2 = C\ wi, A k w k+ \ + B k w k -i = C k w k k G {2, . . . , n - 1}, B n w n _i = C n w n (C.5) 
C\ = —C\ + e. Finally solving the recursion relations we obtain compact expressions for the factors, i.e. 

m— 1 Af l—l m— 1 

n x+eh* 



w. 



i=i 



k=l 



E 

ife>ifc-i>---->ji=2, jk-jk- 



C\ . . . Cj^B^Cj^i . . . Cj k - 2 B jk Cj k+ i ...C, 



m— 1 



m = 21 or m = 21 — 1. 



iM Ai . . . Aj^A^A^+i . . . Aj k _ 2 Aj k Aj k+ i . . . A m „i 

(C.6) 
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(B) It is also worth discussing briefly the cyclic representation. In this case the structure of the 
boundary operator becomes more involved and its diagonalization a more intriguing task. In particular, 
the boundary operator is now of the form: 



p-i 



Qi = ^2 ^ k ek k+i + e v i + X/ Bk ek k ~ i + Bi ei p + X/ ^ k ekk ~ 



k=i 



k=2 



Similarly we define: 



r . s k r .s-\-k r . s-\-k r .s h 

Ak — q 2 : — , r>k — <1 2 _ _ ± , C k —xiq 



q-q- 



and the eigenvalue problem in a matrix form is then written as 
/ Ci Ai Bi \ 



B2 C2 A2 

Bk Ck A k 



Bji—i C n —i A n -i 

B n Cn J 



( W\ \ 

\ w n j 



= e 



( Wl \ 

W2 

\ W n J 



(C.7) 



(C.8) 



(C.9) 



and finally to identify the spectrum and the corresponding eigenstates one has to solve the subsequent 
set of recursion relations: 



A k w k+ i + B k w k -i = C k w k , ke {!,..., p}. 



(CIO) 
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